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Based on a sample of 1.31 x 10® J/'i/' events collected with the BESIII detector at the BEPCII 
collider, Dalitz plot analyses of selected 79,625 rj —>■ 7r"'"7r“7r® events, 33,908 rj —>■ 7r®7r®7r® events and 
1,888 T]' ^ 7r®7r®7r® events are performed. The measured matrix elements of r; —are in 
reasonable agreement with previous measurements. The Dalitz plot slope parameters of r; —>■ 7r®7r®7r® 
and rj' 7r®7r®7r® are determined to be —0.055±0.014±0.004 and —0.640±0.046±0.047, respectively, 
where the first uncertainties are statistical and the second systematic. Both values are consistent 
with previous measurements, while the precision of the latter one is improved by a factor of three. 
Final state interactions are found to have an important role in those decays. 

PACS numbers: 13.66.Be, 14.40.Be 


I. INTRODUCTION 

Since the electromagnetic contribution to the 
isospin violating decays rj/rj' —>■ Stt is strongly sup¬ 
pressed the decays are induced dominantly 

by the strong interaction. Therefore, they offer a 
unique opportunity to investigate fundamental sym¬ 
metries and measure the u — d quark mass differ¬ 
ence. At the tree level of chiral perturbation the¬ 
ory (ChPT), the predicted decay width of ry —?► 
7 r+ 7 r“ 7 r° is about 70 eV, which is much lower than 
the experimental value of 300± 11 eV To explain 
this discrepancy, considerable theoretical effort has 
been made, including a dispersive ajmroach Q and 
non-relativistic effective field theory Q . Recently, it 
was found that higher order terms in ChPT at next- 
to leading order (NLO) Q and next-next-to lead¬ 
ing order (NNLO) @ are crucial for a comparison 
with experimental results, where tttt re-scattering 
between the final state pions is present. 

To distinguish between the different theoretical 
approaches, precise measurements of the matrix el¬ 
ements for rj —>■ TT+TT^TT® and the decay width are 
important. For the three-body decay r] —>■ 
the decay amplitude square can be parameterized 

as[ig 

\A{X, r)p = N{1 + aY + bY"^ -b cX -b dX'^ 
+eXY + fY^ + ...), 


factor. The coefficients a, 6 , c,... are the Dalitz plot 
parameters, which are used to test theoretical pre¬ 
dictions and fundamental symmetries. For example, 
a non-zero value for the odd powers of X, c and e, 
implies the violation of charge conjugation. 

The Dalitz plot distribution of 77 —>■ 7 r+ 7 r“ 7 r° has 
been analyzed previously by various experiments @ . 
Using a data sample corresponding to about 5 x 10® 
r] mesons produced in e“'"e“ —> (j) —> 777 reac¬ 
tions, KLOE provided the most precise measure¬ 
ment, where the Dalitz plot parameters c and e are 
found to be consistent with zero within uncertain¬ 
ties, and / was measured for the first time. Most re¬ 
cently, the WASA-at-COSY collaboration analyzed 
77 —>■ TT+TT^TT® based on a data sample correspond¬ 
ing to 1.2 X 10^ 77 mesons produced in pd —J-^He 77 
reactions at 1 GeV m- The results are in agree¬ 
ment with those from KLOE within two standard 
deviations. 

For 77/77' —>■ 7r®7r®7r®, the density distribution of 
the Dalitz plot has threefold symmetry due to the 
three identical particles in the final state. Hence, 
the density distribution can be parameterized using 
polar variables [l^ 

Z = + (3) 

and the expansion 


where X and Y are the two independent Dalitz plot 
variables defined as 


X = ^{T.+ -T^-), 


( 2 ) 


where Tt^ denotes the kinetic energy of a given pion 
in the 77 rest frame, Q = — m^+ — 777 ^- — m^o is 

the excess energy of the reaction, are the nom¬ 
inal masses from PDG and X is a normalization 


\A{Z)\^ = N{1 +2aZ+ ...), (4) 

where a is the slope parameter, Q = — 3m^o, 

Ti denotes the kinetic energies of each tt® in the 77/77' 
rest frame and X is a normalization factor. A non¬ 
zero a indicates final-state interactions. 

The world averaged value of the Dalitz plot slope 
parameter a = —0.0315 ± 0.0015 Q for 77 — >■ 7r®7r®7r® 
is dominated by the measurements of the Cry stal 
Ball [13, WASA-at-COSY [H and KLOE [I| ex¬ 
periments. Interestingly, the predicted value for a 
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in NLO and NNLO ChPT |,[Tl[i3 is positive, al¬ 
though the theoretical uncertainties are quite large. 

The decay 77 ' —>■ 7 r° 7 r° 7 r° has been explored with 
very limited statistics only. The GAMS-2000 ex¬ 
periment reported the first observation of r]' —^ 
[13 and measured the Dalitz plot slope with 
62 reconstructed events. This result was later up¬ 
dated to be a = —0.59±0.18 [3 with 235 events. In 
2012, the same decay was investigated by BESIII M 
using a data sample of 225 x 10® J/ij} events. The 
branching fraction was measured to be about twice 
as large as the previous measurements, but the 
Dalitz plot slope parameter was not measured. 

In this paper, the matrix elements for 77 —> 
7 r+ 7 r“ 7 r® and 77 / 77 ' —> are measured, where 

the Dalitz plot slope parameter of 77 ' —>■ is 

determined with higher precision than the existing 
measurements. This analysis is performed using a 
sample of 1.31 x 10® J/ip events accumulated with 
the BESIII detector. Radiative J/t/j 777 ^'^ decays 
are exploited to access the 77 and 77 ' mesons. 

II. DETECTOR AND MONTE CARLO 
SIMULATION 

BEPCII is a double-ring e+e" collider working at 
center-of-mass energies from 2.0 to 4.6 GeV. The 
BESIII [13 detector at BEPCII collider, with a 
geometrical acceptance of 93% of 47 r stereo angle, 
operates in a 1.0 T (0.9 T in 2012, when about 
83% of the data sample were collected) magnetic 
field provided by a superconducting solenoid mag¬ 
net. The detector is composed of a helium-based 
drift chamber (MDC), a plastic-scintillator time- 
of-flight (TOE) system, a CsI(Tl) electromagnetic 
calorimeter (EMC) and a multi-layer resistive plate 
counter system (MUG). The charged-particle mo¬ 
mentum resolution at 1.0 GeV/c is 0.5%, and the 
specific energy loss (dE/dx) resolntion is better than 
6 %. The spatial resolution of the MDC is better 
than 130 iim. The time resolution of the TOE is 
80 ps in the barrel and 110 ps in the endcaps. The 
energy resolution of the EMC at 1.0 GeV/c is 2.5% 
(5%) in the barrel (endcaps), and the position res¬ 
olution is better than 6 mm (9 mm) in the barrel 
(endcaps). The position resolution in the MUG is 
better than 2 cm. 

Monte Carlo (MG) simulations are used to esti¬ 
mate backgrounds and determine the detection effi¬ 
ciencies. The GEANT4-based [2l| simulation soft¬ 
ware BOOST [23 includes the geometric and ma¬ 
terial description of the BESIII detector, detector 
response, and digitization models, as well as the 
tracking of the detector running conditions and per¬ 


formance. The production of the J/t/j resonance is 
simulated with KKMC [ 2 ^, [^. w hile the decays 
are generated with EVTGEN for known de¬ 
cay modes with branching fractions being set to the 
world average values Q and by LUNDCHARM 
for the remaining unknown decays. We use a sam¬ 
ple of 1.2 X 10® simulated J/tp events where the J/t/j 
decays generically (‘inclusive MC sample’) to iden¬ 
tify background contributions. The analysis is per¬ 
formed in the framework of the BESIII offline soft¬ 
ware system (BOSS) [13 which takes care of the 
detector calibration, event reconstruction, and data 
storage. 

III. MEASUREMENT OF THE MATRIX 

ELEMENTS FOR THE DECAY 77 ^ 77+77-77° 

For the reconstruction of J/ip —^ 777 with 77 —>■ 
77+77-77® and 77® —^ 77, events consistent with the 
topology 77 + 77-777 s-re selected and the following 
criteria are applied. For each candidate event, we 
require that two charged tracks are reconstructed 
in the MDC and the polar angles of the tracks sat¬ 
isfy |cos0| < 0.93. The tracks are required to pass 
the interaction point within ±10 cm along the beam 
direction and within 1 cm in the plane perpendic¬ 
ular to the beam. Photon candidates are recon¬ 
structed using clusters of energy deposited in the 
EMC. The energy deposited in nearby TOE coun¬ 
ters is included in EMC measurements to improve 
the reconstruction efficiency and the energy reso¬ 
lution. Photon candidates are required to have a 
deposited energy larger than 25 MeV in the barrel 
region (|cos0| < 0.80) and 50 MeV in the endcap 
region (0.86 < | cos0| < 0.92). In order to eliminate 
clusters associated with charged tracks, the angle 
between the directions of any charged track and the 
photon candidate must be larger than 10°. Require¬ 
ments of EMC cluster timing with respect to the 
event start time are used to suppress electronic noise 
and energy deposits unrelated to the event. Events 
with exactly two charged tracks of opposite charge 
and at least three photon candidates that satisfy the 
above requirements are retained for further analysis. 

The photon candidate with the largest energy in 
the event is regarded as the radiative photon orig¬ 
inating from the J/ip decays. For each 77 + 77-777 
combination, a six constraints ( 6 C) kinematic fit 
is performed. The fit enforces energy-momentum 
conservation, and the invariant masses of 77 and 
77+77-77® are constrained to the nominal 77® and 77 
mass, respectively. Events with a from the 6 C- 
kinematic fit 1 ®®® than 80 are accepted for 

further analysis. If there are more than three pho- 


5 



M(7i+7i'7iO) (GeV/c^) 

FIG. 1. Invariant mass spectrum of obtained 

after the kinematic fit without the rj mass constraint ap¬ 
plied. The dots with error bars are for data and the 
shaded histogram is for background events estimated 
from the inclusive MC sample. 


ton candidates in an event, only the combination 
with the smallest Xec retained. To reject possible 
backgrounds with two or four photons in the final 
state, kinematic hts are also performed with four 
constraints enforcing energy-momentum conserva¬ 
tion under the J/^p —>■ 7 r+ 7 r “777 signal hypothesis as 
well as the J/'i/j —>■ 7 r+ 7 r “7777 and J /'0 —>■ 7 r+ 7 r “77 
background hypotheses. Events with a x\c value for 
the signal hypothesis greater than that of the xlc 
for any background hypothesis are discarded. 

After applying the selection criteria described 
above, 79,625 77 —>■ 7 r+ 7 r“ 7 r° candidate events are 
selected. To estimate the background contribution 
under the 77 peak, we perform an alternative selec¬ 
tion, where the 77 mass constraint in the kinematic 
ht is removed. The resulting invariant mass spec¬ 
trum of 7 r+ 7 r“ 7 r°, M( 7 r+ 7 r“ 7 r°), is shown in Fig. [TJ 
A significant 77 signal is observed with a low back¬ 
ground level. The background contamination is es¬ 
timated to be 0 . 2 % from 77 sideband regions, de¬ 
fined as 0.49 < M( 7 r+ 7 r“ 7 r°) < 0.51 GeV/c^ and 
0.59 < M( 7 r+ 7 r“ 7 r°) < 0.61 GeV/c^, in the data 
sample. In addition, a sample of 1.2 x 10® inclu¬ 
sive MG J/t/; decays is used to investigate potential 
backgrounds. Using the same selection criteria, the 
distribution of M( 7 r+ 7 r“ 7 r®) for this sample is de¬ 
picted as the shaded histogram in Fig. [T] No peak¬ 
ing background remains around the 77 signal region. 
The background contamination is estimated to be 
about 0 . 1 %, and is therefore not considered in the 
extraction of the Dalitz plot parameters. 

The Dalitz plot in the variables X and Y is shown 
in Fig. [2] for the selected events. The X and Y 
projections are shown in Fig. [S] For comparison. 
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X 

FIG. 2. Dalitz plot for rj —>■ in the data sample. 

the corresponding distributions obtained from MG 
events with phase space distributed 77 —> 
decays are also shown. The phase space MC distri¬ 
butions of X and Y differ visibly from those in the 
data sample, which indicates there could be large 
contributions from higher order terms in ChPT. 

In order to investigate the dynamics of 77 —>■ 
7 r+ 7 r“ 7 r®, the Dalitz plot matrix elements of the de¬ 
cay amplitude given in Eq. o are obtained from an 
unbinned maximum likelihood fit to the data. To 
account for the resolution and detection efficiency, 
the amplitude is convoluted with a function a{X, Y) 
parameterizing the resolution, and multiplied by a 
function e{X, Y) parameterizing the detection effi¬ 
ciency. Both functions are derived from MC simu¬ 
lations. The sum of two Gaussian functions is used 
for a{X,Y), while e{X,Y) is a quadratic function. 
After normalization, one derives the probability den¬ 
sity function T’(X,Y)^ which is applied in the fit: 

V{X,Y) = 

{\A{X,Y)\^^a{X,Y))e{X,Y) (5) 

{\A{X, r)|2 0 a(X, Y))e{X, Y)dXdY’ 

where A{X,Y) is the decay amplitude of 77 —>• 
7 r+ 7 r“ 7 r® and the integral taken over the Dalitz plot 
(DP) accounts for normalization. 

For the fit, the negative log-likelihood value 

-^event 

-liiC = - (6) 

i=l 

is minimized, where 'P{Xi,Yi) is evaluated for each 
event i, and the sum includes all accepted events. 

We perform two hts to the data. For the hrst ht, 
we assume charge conjugation invariance and we ht 
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FIG. 3. Projections of the Dalitz plot (a) X and (b) Y for 77 —>■ 7r’*’7r“7r'^ obtained from data (dots with error bars) 
and phase space distributed MC events (dashed line). The result of the fit described in the text (solid line) is also 
plotted. 


the parameters for the matrix elements a, b, d and / 
only, while c and e are set to zero. For the second fit, 
we include the possibility of charge conjugation vio¬ 
lation and the latter two parameters are also allowed 
to vary in the fit. 

In the case of charge conjugation invariance, the 
ht yields the following parameters (with statistical 
errors only) 


a = 

-1.128 

± 

0.015, 

b = 

0.153 

± 

0.017, 

d = 

0.085 

± 

0.016, 

f = 

0.173 

± 

0.028. 

The corresponding 

correlation 

matrix 


rameters is given by 



/ 

b 

d 




a 

-0.265 

-0.389 

-0.749 

(8) 


b 

1.000 

0.311 

-0.300 ■ 


U- 


1.000 

0.079/- 

1 


The ht projections on X and F, illustrated as the 
solid histograms in Fig. [3J indicate that the Ht can 
describe the data well. The obtained parameters 
are in agreement with previous measurements within 
two standard deviations. 

If the possibility of charge conjugation violation 
is included in the decay amplitude, the Ht to data 
yields the following results (with statistical uncer¬ 
tainties only) 


a = 

-1.128 

± 

0.015, 

b = 

0.153 

± 

0.017, 

c = 

(0.047 

± 

0.851) X 10-2 

d = 

0.085 

± 

0.016, 

e = 

0.017 

± 

0.019, 

/ = 

0.173 

± 

0.028. 


The corresponding correlation matrix of the Ht pa¬ 
rameters is given by 


/ 

b 

c 

d 

e 


a 

-0.265 

-0.003 

-0.388 

0.001 

-0.749 

b 

1.000 

-0.001 

0.311 

0.016 

-0.300 

c 


1.000 

0.003 

-0.592 

0.003 

d 



1.000 

0.016 

0.079 

\e 




1.000 

-0.007/ 


Compared with the Ht results assuming charge- 
parity conservation, the derived parameters a, b, d 
and / are almost unchanged. The parameters c and 
e are consistent with zero within one standard de¬ 
viation, which indicates that there is no signiHcant 


charge-parity violation in decay 77 —^ Com¬ 

paring the two Hts, the signiHcance of charge-parity 
violation is determined to be only 0.65(7. 

The Ht procedure is veriHed with MC events that 
were generated based on the Dalitz plot matrix ele- 
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merits from the fit to the data. Following the same 
reconstruction and fitting procedure as applied to 
the data sample, the extracted values are consistent 
with the input values of the simulation. 


IV. MEASUREMENT OF THE MATRIX 

ELEMENT FOR THE DECAYS rj 

AND 77' ^ 7rV°7r° 

For the reconstruction of J/4’ —>■ I’n/v' with 
77/77' —>■ 7r°7r°7r*’ and 7r° —>■ 77, events containing at 
least seven photon candidates and no charged tracks 
are selected. The selection criteria for photons are 
the same as those described above for 77 —7r+7r“7r°, 
except the requirement of the angle between the pho¬ 
ton candidates and any charged track. Requirements 
of EMC cluster timing with respect to the most en¬ 
ergetic photon are also used. Again, the photon 
with the largest energy in the event is assumed to 
be the radiative photon originating from the J/if) 
decay. From the remaining candidates, pairs of pho¬ 
ton are combined into 7r° —^ 77 candidates which 
are subjected to a kinematic fit, where the invariant 
mass of the photon pair is constrained to the nomi¬ 
nal 7r° mass. The value of this kinematic fit with 
one degree of freedom is required to be less than 25 . 
To suppress the mis-combination, the decay 
angle 0decay, defined as the polar angle of a photon 
in the corresponding 77 rest frame, is required to 
satisfy |cos 0 decay| < 0 . 95 . From the accepted 7r° 
candidates and the corresponding radiative photon, 
'yTT^TT^TT^ combinations are formed. A kinematic fit 
with seven constraints ( 7 C) is performed, enforcing 
energy conservation and constraining the invariant 
mass of 77 pairs to the nominal 7r° mass. If more 
than one combination is found in an event, only the 
one with the smallest is retained. Events with 
Xyc < ’^0 are accepted for further analysis. 

For 77' —>■ 7r°7r°7r°, backgrounds from J/t/; —>• 
ujTT^TT^ are suppressed by vetoing events with 
IM(777°) —777(^1 < 0.05 GeV/c^, where M(77r°) is the 
invariant mass of the 771^ combination closest to the 
nominal ui mass {rrii ^) @ . Peaking backgrounds for 
the process 77' —>■ can arise from J/t/ —>■ 777' 

with 77' —>■ rjir^TT^. To suppress these backgrounds, 
a 7 C kinematic fit under the J/t/ —>■ xrjTr^Tr^ hy¬ 
pothesis is performed. Events for which the 
value obtained for the background hypothesis is less 
than that obtained for the 77r°7r°7r° hypothesis are 
discarded. In addition, events with an invariant 
mass of at least one 77 pair in the mass window 
\M{n) ~ ™r;| < 0.03 GeV/c^ are rejected. 

For 77 —^ TT^TT^TT^, the invariant mass spectrum 
of is shown in Fig. IH^a). A very clean 77 


signal is observed. The invariant mass spectrum 
of 7r°7r°7r° obtained from the inclusive MG sample 
is also shown, indicating a very low background 
level of 0 . 3 % under the 77 signal. The background 
is also estimated from the data using 77 sideband 
regions ( 0.49 < M{n°Tr°Tr°) < 0.51 GeV/c^ and 
0.59 < M{TT°TT°n°) < 0.61 GeV/c^), and is found 
to be about 1 %. For the determination of the slope 
parameter a, the backgrounds are neglected. 

To improve the energy resolution of the can¬ 
didates and thus the resolution of the Dalitz plot 
variable Z, the kinematic fit as described above 
is repeated with the additional constraint that the 
7r°7r°7r° invariant mass corresponds to the nominal 77 
mass. 

Finally, a clean sample of 33,908 77 —>■ 
events is selected. The distribution of the variable 
Z, defined in Eq. ([ 3 ]), is displayed in Fig. SKb). The 
dotted histogram in the same plot represents the 
MG simulation of phase space events with a = 0 , 
as expected at leading order in ChPT. Due to the 
kinematic boundaries, the interval of 0 < Z < 0 . 7 , 
corresponding to the region of phase space in which 
the Z distribution is flat, is used to extract the slope 
parameter a from the data. 

Analogous to the measurement for 77 —> 
an unbinned maximum likelihood fit is performed 
on the Z distribution of the data to extract the 
slope parameter. The probability density function 
is constructed with Eq. o convoluted with a dou¬ 
ble Gaussian function and multiplied by a first- 
order Ghebychev polynomial to account for the res¬ 
olution cr(Z) and detection efficiency e{Z), respec¬ 
tively. Both the resolution and the efficiency func¬ 
tions are obtained from the phase space distributed 
MG events. The fit yields a = — 0 . 055 ± 0 . 014 , where 
the error is statistical only. In the inset of Fig. IHb) 
the result of the fit is overlaid on the distribution for 
the data. 

For 77' —>■ TT^TT^TT^, the invariant mass spectrum 
of 7r°7r°7r° is shown in Fig. EKa), where an 77' sig¬ 
nal is clearly visible. The analysis of the J/ij) inclu¬ 
sive decay samples shows that the dominant back¬ 
ground contribution is from 77' —>■ rjir^Tr^. Addi¬ 
tional backgrounds are created by J/t/j decays to the 
same final state, e.g., J/t/ with w —?► 777°. 

To evaluate the contribution from 77' —> 

4 X 10 ® J/t/j —>■ 777' events with 77' —>■ 7777 ®tt® are gener¬ 
ated. The 77' decay dynamics are modeled according 
to the results of the Dalitz plot analysis given in 
Ref. [2^. The invariant mass spectrum of 7r®7r®7r® is 
also shown in Fig. EJa), where the number of events 
is scaled to the number of J/t/j events in the data 
sample, taking into account the branching fractions 
of J/t/ —>• 777' and the subsequent decays. Other 



0.45 0.5 0.55 0.6 0.65 

M(7i°7tV) (GeV/c^) 



FIG. 4. (a) Invariant mass spectrum of for rj obtained from data (dots with error bars) and 

estimated from the inclusive MC sample (shaded histogram), (b) Distribution of the kinematic variable Z for rj —>■ 
7 r°7r°7r° obtained from data (dots with error bars) and phase space distributed MC events, where the Z distribution 
is flat from Z = 0 to Z ^ 0.76 and then drops to zero at Z — 1 (dashed line). The inset shows the Z range which is 
used for the fit to extract the slope parameter a. Overlaid on the data is the result of the fit (solid line in the inset). 
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FIG. 5. (a) Invariant mass spectrum of for rj' obtained from the data (dots with error bars), 

estimated from the inclusive MC sample (shaded), and rj' —>■ rj'K^'K^ MC events (hatched), (b) Distribution of the 
kiirematic variable Z for tj' —>■ obtained from the data (dots with error bars), phase space distributed MC 

events (dashed line), rj' sideband regions (shaded) and tj' —>■ rjTv^TT^ MC events (hatched). The result of the fit (solid 
line) is overlaid on the data in the insert. 


background contributions {e.g. from J/ip ^ cutt^tt^) 
are estimated from the data sample using the rj' side¬ 
band regions, defined as 0.845 < M(7r°7r°7r°) < 0.88 
GeV/c2 and 1.008 < M{7r°7r°Tr°) < 1.043 GeV/c^ 
(Fig. EKa)). The total background contamination is 
estimated to be 11.2% in the rj' signal mass region 
(0.92 < M(7r°7r°7r°) < 0.99 GeV/c^). 

After requiring the invariant mass of 7r°7r°7r° to 
be in the rj' signal mass region, the distribution of Z 
is shown in Fig. EKb). The MC simulation of phase 
space events clearly deviates from the data. Analo¬ 


gous to 77 —?► 7 r° 7 r° 7 r°, the slope parameter a is deter¬ 
mined from an unbinned maximum likelihood fit to 
data in the range 0 < Z < 0.45 with 1,888 events, 
taking into account the detection efficiency and reso¬ 
lution. The background estimated from rj' —> 

MC events and the ij' sideband regions is accounted 
for by subtracting the likelihood for these events 
from the likelihood for data. The normalization of 
background contribution is fixed at its expected in¬ 
tensity. 

The fit yields a slope parameter a = —0.640 ± 

































0.046, where the error is statistical only. The result 
of the fit is overlaid on the Z distribution for the 
data in the inset of Fig. [5l[b). 

V. SYSTEMATIC UNCERTAINTIES 

Various sources of systematic uncertainties on the 
measured Dalitz plot matrix elements have been in¬ 
vestigated. These include uncertainties due to the 
efhciency parameterization and uncertainties aris¬ 
ing from differences in the tracking and tt^ recon¬ 
struction between the data and MC samples. For 
the measurement of a for ij/rj' —>■ addi¬ 

tional uncertainties due to the fit range and mis- 
combination are considered. Uncertainties for a. due 
to the background estimation for r\' —^ are 

also assigned. All the above contributions are sum¬ 
marized in Table m where the total systematic uncer¬ 
tainty is given by the quadratic sum of the individ¬ 
ual errors, assuming all sources to be independent. 
Assuming the correlation factor between each sys¬ 
tematic errors is 1 , then the correlation matrix for 
systematic errors of r/ —)■ is 


/ 

b 

d 


a 

-0.71 

0.99 

-0.97 

b 

1.00 

-0.73 

0.54 

\d 


1.00 

-omJ 


In the following, the estimation of the individual un¬ 
certainties are discussed in detail. 

To estimate the uncertainty due to efficiency 
parameterizations, we perform alternative fits by 
changing the description of the efficiency from poly¬ 
nomial functions to the average efficiencies of local 
bins. The change in the obtained values for the ma¬ 
trix elements from the alternative fits with respect 
to the default values is assigned as the systematic 
uncertainty due to the efficiency parameterization. 

Differences between the data and MC samples for 
the tracking efficiency of charged pions are investi¬ 
gated using J/'0 —>• ppTT~^'K~ decays. A momentum- 
dependent correction is obtained for charged pi¬ 
ons reconstructed from MC events. Similarly, a 
momentum-dependent correction for the effi¬ 
ciency in the MC sample is obtained from J/if) ^ 
7 r+ 7 r“ 7 r° decays. The fits to extract the matrix ele¬ 
ments are repeated as described above, taking into 
account the efficiency correction for charged pions 
and TT^. The change of the matrix elements with 
respect to the default fit result is assigned as a sys¬ 
tematic uncertainty. 

The slope parameter a for 77 / 77 ' —is ex¬ 
tracted from a fit to the data in the kinematic re¬ 
gion where the Z distribution of phase space is flat. 
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By altering the fit range to 0 < Z < 0.65(0.68) 
for ?7 —>■ 7 r° 7 r° 7 r° and 0 < Z < 0.43(0.45) for 
77 ' —>■ 7 r° 7 r° 7 r° and repeating the ht to the data, the 
larger changes in a with respect to the default fits 
are noted and assigned as the systematic uncertain¬ 
ties. 

Mis-reconstruction of candidates in true signal 
events can lead to a wrongly reconstructed position 
of the event on the Dalitz plot, and therefore affect 
the fitted parameters. Using signal MC, the possi¬ 
ble mis-combination of photons has been studied by 
matching the generated photon pairs to the selected 
7 r° candidates. The fraction of events with a mis- 
combination of photons is 5.4% for 77 —and 
0.95% for 77 ' —>• 7 r° 7 r° 7 r°, respectively. Applying the 
fit to the truth-matched simulated events only, the 
impact on the fit parameters is found to be 2 . 8 % for 
77 —>■ 7 r° 7 r° 7 r° and 1 . 0 % for 77 ' —>■ 7 r° 7 r° 7 r°, respectively. 
This is taken as the systematic uncertainty. 

In the determination of a for 77 ' —>■ back¬ 

ground contributions are estimated from MC simu¬ 
lations for 77 ' —>■ 7777 ^ 77 ° and 77 ' sideband regions. For 
the peaking background from 77 ' —>■ 7777 ^ 77 ®, the un¬ 
certainties of the branching fractions for J/t/ —>■ 777 ' 
and 77 ' —>■ 7777 ^ 77 ° taken from Ref. @ are consid¬ 
ered. In addition, an alternative set of matrix ele¬ 
ment parameters for 77 ' 7777 ^ 77 ° as reported by the 

GAMS-47r collaboration in Ref. [ 2 ^ is used in the 
MC simulation. The uncertainty from non-peaking 
backgrounds is estimated by varying the sideband 
regions to 0.723 < M(7r°7r°7r°) < 0.758 GeV/c^ and 
1.063 < Af ( 77077077 °) < 1.098 GeV/c2. 

In order to estimate the impact from the differ¬ 
ent resolution of Dalitz plot variables between data 
and MC sample, we perform alternative fits in which 
the resolution is varied by ± 10 % and find that the 
change of the results is negligible, as expected. 

VI. SUMMARY 

Using 1.31 X 10° J/ 7/7 events collected with the 
BESIII detector, the Dalitz plots of 77 —>■ tt+tt^ttO 
and 77 / 77 ' —)■ TT^TT^TT^ are analyzed and the corre¬ 
sponding matrix elements are extracted. 

In the case of charge conjugation invariance, the 
Dalitz plot matrix elements for 77 tt+tt^ttO are 
determined to be 

a = -1.128 ±0.015 ±0.008, 

b = 0.153 ±0.017 ±0.004, 

d = 0.085 ±0.016 ±0.009, 

/ = 0.173 ± 0.028 ±0.021, 

where the first errors are statistical and the second 
ones systematic, here and in the following. In Fig. [5] 
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TABLE I. Summary of systematic uncertainties for the measurements of the matrix elements (all values are given in 
%). 


Source 

a 

b 

d 

/ 

0(77 — >■ 7r°7r°7r°) 

a(r]' 7r°7r°7r°) 

Efficiency parameterization 

0.6 

1.7 

10.4 

11.7 

0.4 

0.1 

Tracking efficiency 

0.1 

0.6 

0.2 

0.2 

- 

- 

7r° efficiency 

0.1 

2.0 

1.6 

1.3 

3.7 

1.6 

Fit range 

- 

- 

- 

- 

3.7 

3.4 

7r° mis-combination 

- 

- 

- 

- 

2.8 

1.0 

Background subtraction 

- 

- 

- 

- 

- 

6.2 

Total 

0.7 

2.7 

10.5 

11.8 

6.1 

7.3 


our measurement is compared to previous measure¬ 
ments and theoretical predictions. Our results are in 
agreement with the two most recent measurements, 
and consistent with the predictions of the dispersive 
approach and ChPT at NNLO level. 

To investigate the charge conjugation violation in 
77 —^ TT+TT^TT*^, the matrix elements c and e have been 
determined from a fit to the data. The obtained 
values are consistent with zero, while the other pa¬ 
rameters are found to be consistent with those ob¬ 
tained from the fit assuming charge conjugation in¬ 
variance. No significant charge symmetry breaking 
is observed. 

After taking into account the systematic uncer¬ 
tainties, the slope parameter a for 77 —^ 7 r*^ 7 r*^ 7 r'^ is 
measured to be —0.055 ± 0.014 ± 0.004. A compar¬ 
ison to previous works, illustrated in Fig. [7Ka), in¬ 
dicates that the BESIII result is compatible with 
the recent results from other experiments and in 
agreement with the prediction from ChPT at NNLO 
within two standard deviations of the theoretical un¬ 
certainties. 

The Dalitz plot slope parameter for 77' —>■ 7r°7r°7r° 
is measured to be a = —0.640± 0.046± 0.047, which 
is consistent with but more precise than previous 
measurements (Fig. EKb)). The value deviates sig¬ 
nificantly from zero. This implies that final state 
interactions play an important role in the decay. Up 
to now, there are just a few predictions about the 
slope parameter of 77' —?► In Ref. [ 2 ^, the 

slope parameter is predicted to be less than 0.03, 
which is excluded by our measurement. More re¬ 
cently, using a chiral unitary approach, an expan¬ 
sion of the decay amplitude up to the fifth and sixth 
order of X and Y has been used to parameterize 
the Dalitz plot of 77' —^ [s^. The coefficient, 

which corresponds to a in this paper, is found to be 
in the range between —2.7 and 0.1, consistent with 


our measurement. 
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. * Theoretical predictions without error. **BSE denotes Bethe-Salpeter Equation. 
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